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!>• ! Abstract 

. 

t—i ' In this letter, dark energy is obtained using dual roles of the Ricci 

scalar (physical as well as geometrical role of the Ricci scalar curvature) 

c| . in a natural way without any input in the theory for it. In contrast 

Or 

^ | to other models, where an idea for dark energy is introduced in the 

5h ' 

^ ■ beginning of the theory, here it emerges from the gravitational sector 

03 ■ 

spontaneously Dark energy density, obtained in this model, decreases 

^ ■ 

with the scale factor a(t) as ~ a(i)~ 3 ( 1+w ( de )) (where W(d e ) > — 1 is 

the equation of state parameter for dark energy ) explaining fall of 

its value from 1.19 x 10 75 GeV 4 at Planck scale to its present value 

suggested by astronomical observations. Interestingly, apart from dark 

energy, dark radiation and dark matter density are also obtained from 

the gravity sector. Moreover, time for transition from deceleration to 

acceleration of the universe is evaluated. It is found that if dark energy 

obeys equation of state for generalized Chaplygin gas and barotropic 

fluid simultaneously, after tansition time, the scale factor will evolve 

as cosh(t — t*)(t and £* being time and transition time respectively), 

l 
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but after a finite time expansion will become de-Siiter like. PACS no. 
98.80.Cq 



f. Astronomical observations, during the last few years, provide com- 
pelling evidences in favor of accelerating universe at present, which is 
caused by dominance of dark energy (DE) [1, 2, 3, 4, 5]. Observation 
of 16 Type la supernovae (SNe la) by Hubble Space Telescope further 
modifies these results and shows evidence for cosmic deceleration pre- 
ceding acceleration in the late universe [6]. So, DE has a significant 
role in cosmic dynamics of the current universe. The simplest candi- 
date for DE is supposed to be the cosmological constant A which is 
very high in the early universe, but there is no mechanism to bring it 
to present value without fine-tuning. Alternatively, to explain its decay 
from a very high value, in the early universe, to its present extremely 
small value, many models [7, 8] were suggested in which A is envisaged 
as a slowly varying function of cosmic time. Apart from dynamical 
A as DE, other models are fluid-dynamics models, where barotropic 
fluid is its source with or without dissipative pressure , Chaplygin 
gas(GC) and generalized Chaplygin gas (GCG)models [9, 10, 11, 12]. In 
field-theoretic models, the most natural ones are models, where DE is 
caused by scalars. These are quintessence models [13], k-essence mod- 
els [14], tachyon models [15] and phantom models [12, 16]. In these 
models, non-gravitational lagrangian density for exotic matter giving 
DE is added to Einstein-Hilbert term in the action. Recently, in a 
different appraoch, non-gravitational term is replaced by gravitational 
term being non-linear in Ricci scalar R, which stems modified gravity 
[17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28]. 
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In all these models, a non-gravitational or gravitational term for 
DE is added in the theory a priori. In this sense, these models are 
phenomenological. Cosmology demands a DE model, where DE density 
P(dc) emerges spontaneously from a basic theory. With this motivation, 
here,p(de) is derived from the gravitational sector without an input for 
it. Here also, a quardatic term of R is added to Einstein-Hilbert term 
in the action, but not as lagrangian density for DE. It is unlike the 
approach in [17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28]. In the pesent 
model, DE emerges spontaneously due to presence of term quardatic 
in R manifesting dual role of the Ricci scalar as a physical field as well 
as geometry. 

Natural units (h — c — 1) are used here with GeV as the fundamental 
unit, where h and c have their usual meaning. In this unit, it is found 
that lGeV" 1 = 6.58 x l(T 25 sec. 

2. The action for higher- derivative gravity is taken as 



where R is the Ricci scalar curvature, G = Mp 2 (M P = 10 19 GeV is 
the Planck mass). Moreover, dimensionless a is a scalar depending 
space-time coordintaes. 

The action (1) yields gravitational field equations 



-2(V„V„a - g^Oa)R - 4(V M aV^ - g^V a aV a R) = (2) 



using the condition 5S g /5g^ u = 0. Here, denotes covariant derivative 
and the operator □ is given as 




(1) 



1 



(R„v - \g^R) - a(2V^ u R - Ig^UR - \g^ v R 2 + 2RR^) 



167tG 




(3) 
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with 11,1/ — 0, 1, 2, 3 and as metric tensor components. 

Taking trace of eqs.(2), it is obtained that 

2 

DR + m 2 R + -S/ v aS/ u R = 0, (4a) 
a 

where 

m 2 = — — H (46) 

967rGa(t) a y ' 

with a 7^ to avoid the ghost problem. Here, overdot gives derivative 
with respect to time t. 

Eq.(4a) shows that the Ricci scalar R behaves as a physical field 
also with (mass) 2 depending on G, in addition to its usual role as a 
geometrical field [29, 30, 31]. 

Experimental evidences support spatially homogeneous flat model of 
the universe [32]. So, the line-element, giving geometry of the universe, 
is taken as 

dS 2 = dt 2 - a 2 (t) [dx 2 + dy 2 + dz 2 } (5) 

with a(t) as the scale factor. 

In the space-time, given by eq.(5), eq.(4a) is obtained as 

- /a a\ ■ r 1 a dd^i 

R+(3- + 2-)R-\ + - + 3— \R = 0, 6 

V a a/ (-9D7rGo;(t) a a a J 

In most of the situations, for example, radiation model, matter- 
dominated model, and accelerated models, we have a(t) as a power- 
law solution yielding R as the power-law function of a(t) . So, it is 
reasonable to take R as 

R = 4- (7) 

R n V I 

with n^O being a real number. Moreover, a is taken as 



a = Da r , 



(8) 
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where r is a non-zero real number. Advantage for taking a in this form 
is mentioned below. Here A and D are constants. 

R and a, given by eqs.(7) and (8) respectively, satisfy eq.(6), if 

_ + (2 + r-n - =-——— — , 9 

a \a/ 96nGD(n — r)a r 



which integrates to 
-a\2 C 



96%GD(n -r)a r ' 

(10) 



.a) ~ a 2 I 3 -( n - r )] 48%GD(n - r)[(6 - r) - 2(ra - r)]a r 
with C being an integration constant. Eq.(10) gives dynamics of the 

universe, which is the Friedmann equation 

'd\2 8nG 



with 



3C 

Pdc ~ 8nGa^-("-r )1 ( } 



and 



Pr(m) 



8nG 



48nGD(n - r)[(6 - r) - 2(n - r)]a r J ' ^ 12 ^ 
It is interesting to note that a(t) oc a r , gives energy density in the 

known form oc a~ r . 

p r (m), given by eq.(12b), takes the known form of energy density in 
two cases (i) r = 4 and (ii) r = 3. When r = 4, it takes the form 
of radiation density p r ~ a~ 4 . It is in the form of matter (pressureless 
fluid) density p m ~ a~ 3 , when r = 3. 

Pdc, given by (12a), emerges from the gravitational sector without 
using a source for it. So, it is recognized as DE density. 

Conservation equation for DE 

P(dc) + 3-p ( de)(l + W(de)) = 0, (13) 
CL 

where equation of state (EOS) parameter W(d e ) = P(de)/P(de) with p^e) 
and p(de) being isotropic pressure and density for DE respectively. 
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Eq.(13) yields 

2[3-(n-r)] =3(l + w (dc) ), (14) 
Thus, eqs.(12a) and (14) imply 

Pde ~ 87rGa^ 1+ ^y (15) 
WMAP data give the present value of DE to be p® de ^ = 0.73p cr ., 

where p cr . = 3H%/8irG with H = 100/ikm/secMpc = 2.33 x l(T 42 GeV 

and h = 0.68 (having maximum likelihood) [34, 35]. Using these values, 

eq.(17) is obtained as 

Pde = 0.73p cr . — . (16) 



a 

3 A/T2 



At Planck scale, p( de ) — s^g^P' so e Q-(24) yields 

(M*'«' = W = s . (6xl0 ... (17) 
\a P J 0.73 ' v ; 

where dp is the scale factor at Planck time tp — Mp 1 . 
2(a). When r = 4, eq.(12b) is obtained as 

Pr = 8^G .72nGD(l - w dc )| - 1 + 3w dc |a 4 ] ' ^ 
which looks like energy density for radiation emerging from gravity sec- 
tor. This type of term also arises in brane-gravity inspired Friedmann 
equation, which is called as 'dark energy density' [33]. So, like brane- 
gravity, here also p r (given by eq.(18)) is termed as 'dark radiation 
density'. 

Connecting eqs.(ll), (12a) and (16), it is obtained that 

(^) 2 = 4 + 0.73^(^) 3,I+ •' < " , (19) 
\a/ a 4 V a / 

with E = l/{72nGD(l - w dc )| - 1 + 3w dc |}. 
For p r > p( de ) , eq.(19) reduces to 

E 
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which integrates to 

a 2 = a 2 (0) + 2t\[E (21) 

with a(0) = ait = 0). It shows deceleration as a < 0. 

Investigations start here from the Planck scale (being the fundamen- 
tal scale). Using ap = a(t = tp) with Planck time tp = Mp 1 in eq.(18), 
C is evaluated as 

C = ^(a 2 P -a 2 (0)) 2 M P ~^a P M P . (22) 
as a(0) < ap. Thus, 

a 2 ~ a P M P i. (23) 
Connecting eqs.(16), (17), (18) and (22), it is obtained that p r > p dc 

if 

^y-3w de ^ 5 ^ ^ iQ _ ll2 ^ ^ ^ ^ 10l21)4 /3(l + w de )_ (24) 

For w dc = —0.89, this inequality yields 

1 + ^ = ^> 10 365 , (25) 
a 

where z is the red-shift. 

It is found that universe decelerates when p r > pde- Moreover, uni- 
verse accelerates when DE dominates. It means that transition from 
dceleration to acceleration is possible when transition from p r > p^ c 
to p r < pdc takes place. According to (25), it is possible at red-shift 
10 365 , which is extremely large. 16 Type Supernova suggest red-shift 
for this transition as z = 0.46±0.13 [6]. Thus this case (for (n — r) = 1) 
contradicts experimental results. So, it is ignored here onwards. 

2(b). When r = 3, eq.(12b) looks like 
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with B = 01g „„ 1 — r, — r . The energy density (26) has the form of 

2167rG(l-w do )|w dc | OJ 

density for pressueless matter emerging from gravity. So, it is termed 
as dark matter density (DMD). The prtesent value of DMD is p® m = 
0.23p cr (critical density p cr is defined above). So, (26) is obtained as 

p m = 0.23p cr (^) 3 . (27) 

Now, p m > pdc when red-shift z is given as 

/73\V3|w de | 

*>(™J -1- (28) 



,23 

Pm < Pdc for 



/73\i/3|w dc | 

*<(™) - 1 - (29) 



,23, 

Inequalities (28) and (29) show that transition from p m > p dc to 

Pm < Pdc at 

Observations yield — 1 < w dc < —0.82 for quintessence DE. Using 
these values of w dc , (30) yields 

0.47 <z*< 0.599, (31) 

which is supported by the range 0.33 < z* < 0.59 given by 16 Type 
Supernova observations. 

Using (16) and (27), Friedmann equation (11) is obtained as 



) 2 = i/ 2 [0.23(^) 3 + 0.73(^) 3(1+Wdc) ] (32) 



'a 
.a 

using definition of p° cr given above. 

When z > z*,p m > p dc , so (32) reduces to 

2 _ /fin \ 3 
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which integrates to 



a(t) = a d \l + 0.72H (^) 3/ \t-t d ) 



1 2 /3 



(34) 



with ad = a(td) being a constant, a(t), given by (34), shows deceleration. 
Further, for z < z*,p m < pdc, so (32) reduces to 



0.73#n 



2 /Oq\ 



3(l+w dc ) 



a J 



which is integrated to 

r + 3(l+w de ) J/ov/ ^ ( aox3q + » J .)/2 (t _ tj , 1 



a(t) = a* 







(35) 



(36) 



2 ' Va*/ 

with a* = a(t*). (36) shows acceleration. Here t* is the time for 

transition from deceleration to acceleration. Using a = a(t ) in (36), 
we obtain that 



tn — t* 



1 + 3(1 + Wdc) jW073 



- ; - /a*\ 3 ( 1+Wd «V 2 
1- - 
Va / 



(37) 



Connecting (30) and (37), t* is evaluated as 

U = 0.627t = 8.59Gyr (38) 
for w dc = —0.89 and present age of the universe t = 13.7Gyr [34]. 

3. So far, we have considered DE (obtained here) as a barotropic fluid 
obeying EOS 

Pdc = w dc p de . (39a) 

Due to negative pressure and being only fluid with supersymmetric 
generalization , GC and GCG have been strong candidates for DE. 
But experimental results support GCG [10, 11]. So, GCG model is 
preferred. 

In what follows, consequences of DE fluid (obtained here) are ex- 
plored if it behaves as GCG [9, 10, 11, 12] and barotropic fluid simulta- 
neously. This type of situation is explored in [12] for non-gravitational 
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phantom DE. The equation of state (EOS) for GCG is given as 



Pde = 



13 
Pde 



(396) 



where < (3 < 1 for GCG , (5 — 1 for CG and M is a constant. 
(39a) and (39b) yield 



w dc = - 1 — 1 



/ M y+p 



(40) 



This equation shows dependence of Wd c on pd e varying with a(t). Now, 
using eq.(16) for pd c with variable Wd e , conservation equation (13) looks 
like 



3w, 



de 



3 /d 
Ina \a 



using eq.(40) for w de . Here, 

B = 0.73 Pcr af +Wdc) . 

Eq.(41a) integrates to 

" fl 3(l+/3) + Ui ' 

where C is an integration constant. 
So, eqs.(16), (41b) and (42) yield 

l/(l+/3) 



)lna 

(41a) 
(416) 

(42) 



Pde = B 



_ a 3(l+/3) ^ \ B ) 

Connecting eqs.(40) and (43), it is obtained that 

vb) 



w dc 



(#) 



1 a 3(l+/3) + 1 ~ I 



(43) 



(44) 



At a = a ,w dc = w° c , so eq.(44) yields 



/My+p 
\~b) 



Connecting eqs.(44) and (45a), we obtain 

l+w° \ /a \3(i+/3) 



w dc 



w 



de 



a) 



- 1 



(456) 
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Eqs.(43) and (45a) yield 

BCW 1 *® r / a \ 3(i+/3) / w° dc \ i ' ! - ' ' 
a) \l+w° de . 



11 



Pdc 



(46a) 



At a = a , p dc = P° c = 0.73p cr , so 



p dc = 0.73p c 



o ,/a \W) 



w. 



de 



V(l+/3) 



(466) 



setting C = a- 3w S=( 1+ ^. 

The constant M, in (40), is evaluated as 

M 1+ ? = -< c (p° dc ) 1+ P 



(47) 



using w° e = Wde(to)- So, connecting (40), (46b) and (47), it is obtained 
that 



w de = w de 



s.)(t) 



a \3(i+/3) 



w 



dc 



(48) 



In the changed circumstances, p m < p dc when 



23 /oq\3 
73 V a) 



/a \W) 
l+w H J — ) -w, 



dc 



l/(l+/3) 



(49) 



This inequality is obtained using (27) and (46b). It shows that, for 
GCG , transition from p m > p dc to p m < p dc may take place when 



/a \W) 

— = ~ w d 

Va* / gcg 

Taking w° de = -0.89 in (50), 



^*(GCG) 



(73) - (1+W ° J 



^] - 1 = 0.64 



(50) 



(51) 



for P = 0.04. This value of red-shift for transition is very near to 
upper limit 0.59 of z, given in [6]. So, in GCG case, p m < p dc when 
z < z* — 0.64. Now, using (46b) and definition of p cr , Friedmann 
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Equation looks like 

(-) = 0.73^[(l + wS o )(f) -w»«] 

(1 + 2/3) / (l+w° e ) ^ ) 
4(1 + /3)H KJ J a ° 



(52) 



which is integrated to 

-y/T+2PcoshO - 1] i/3(i+/») 



a(t) = a* (G CG) [ — -j===— - — J , (53a) 

where 

9 = 3v / 0^3i/o(l +/?)|w° c | 1/2(1+/3) (^-^(gcg)). (536) 
Using f3 = 0.04 and = 0.64, £*(gcg) is evaluated as 

^(gcg) - 0.8t = 10.9Gyr. (54) 
Connecting (53a) and (53b), it is obtained that 

* = yfimH^V*™ \ f^ tmM J ■ (55a) 
a La/1 + 2p — sech0\ 

It shows that 

- ~ v / 0J3// |w° e | 1 /2(i+/?) (556) 

(X 

for = 12.2. Using from (53b), The corresponding time t& s , for 
6 = 12.2, is obtained as 

U s = U GCG) + 4 -^^\< c \- 1/2(1+ " ) - (56) 
Using (3 = 0.04 and w° c = —0.89 in (50), t ds is evaluated as 

Us = t*(GCG) + 4.88t ^ 5.68t = 77.78Gyr, (57) 
where £*(gcg) is given by (54). 
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For t > ids, tanhO = 1 and coshO ~ y , so a(t) (given by (53a)) looks 



like 

v'l -2 iroshO - 2-|V : HH n 



a(t) 



► (GCG) 



VI + 2/3 -2 
.^IW 1 1/3(1 ; 



Fe V0l3H 1 w» c | V2(i+« (*-t. (GOG) ) 

(58) 



with 6 1 , given by (53a). Moreover, 







= 0J3if o Vde| 1/2(1+/?) 

GCG U 



where F = a* (G cG) 
for t>t ds , 

Pde(GCG) -O.TSpcrlw^l 1 ^ 1 ^. (59) 

(58) shows that ,in GCG case,when t > tds, expansion becomes de- 
Sitter like and and DE density acquires a constant value slightly less 
than present DE density. In non-GCG case, p dc — > as t — > oo. From 
(48), it is found that —0.63 > w dc > —1 taking (3 and w° e given above. 
But for t > t ds ,w dc ~ —1. 
Moreover, (53a) yields 

sechO i - 1 . . 

VTT2^J 1 ; 

This equation shows that as t increases from t = i*(GCG), ) de- 
creases and becomes ~ 0. 73iJ,f | w° c 1 1 / 2 ( 1 +^ ? ) for t > t ds . It means that, in 
GCG case, expansion is faster around transition time, but slows down 

In non-GCG case, a(t) is given by (36), which yields 

(*) =0.73^[l- 3(1 + Wdc) lf^ 3(1+Wdc) . (61) 

W non-GCG °L 2 JVa/ 

Thus, it is found that acceleration increases with growing scale factor 
a(t) in both GCG and non-GCG cases. But, in GCG case, a oc a 
and a oc a 1-3 ( 1+Wde ) f or non-GCG case. It means that, in GCG case, 
acceleration grows faster with increasing a(t) compared to the non- 
GCG case if w de > — 1. Its reason is explained in the following way. 



14 S.K.SRIVASTAVA 

Accelerated expansion is caused by DE, characterized by negative 
pressure violating the strong energy condition, which gives cosmic push 
due to its reverse gravity effect. In non-GCG case, — p decreases with 
increasing scale factor (as — p oc p), but it increases with a(t) for GCG 
case and attains a constant value in finite time. So, GCG dark energy 
gives more cosmic push causing more acceleration compared to non- 
GCG dark energy, 

4. In this letter, it is found dark energy emerges spontaneously from 
gravity sector using dual nature of the Ricci scalar as a geometrical field 
as well as physical field (mentioned above). Quintessence dark energy 
density, obtained here, falls from very high value 1.19 x 10 75 GeV 4 at 
the Planck scale to its current value 0.73p cr . ~ 2.18 x 10~ 47 GeV 4 . At 
the transition time i*, which is 8.59Gyr for w^c = —0.89, dark energy 
begins to dominate giving a sudden jerk to the universe. As a result, 
a transition from deceleration to acceleration takes place showing the 
reverse gravity effect of dark energy. It is also found that dark radia- 
tion is washed away in the late universe. Transition from dominance 
of DM to dominance of DE is obtained at red-shift 0.47 < < 0.599 
(for —0.82 > Wdc > — 1) causing transition from deceleration to accel- 
eration. Interestingly, z*, obtained here, is consistent with the value 
of given by 16 Type SNe la observations. Universe accelerates after 
t > t* with power-law expansion with vanishing energy density, when 
t — > oo. But, if DE fluid behaves as GCG and barotropic fluid si- 
multaneously, accelerated expansion is faster due to more cosmic push 
compared to non-GCG case. It is found that, in GCG case, expansion 
becomes de-Sitter like after a finite time with finite energy density. 
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